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We perform a fit to the 2006 B — » tvK data and show that there is a disagreement with the 
standard model. That is, the B — > nK puzzle is still present. In fact, it has gotten worse than in 
earlier years. Assuming that one new-physics (NP) operator dominates, we show that a good fit 
is obtained only when the electroweak penguin amplitude is modified. The NP amplitude must be 
sizeable, with a large weak phase. 
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PACS numbers: 13.25.Hw, 11.30.Er, 12.15Ff, 14.40.Nd 

There are four B — > nK decays - B + — ► tt + K° (desig- 
nated as +0 below), B+ -► tt°K+ (0+), B° d -> ir~ K+ 
(—+) and B d — > tt°K° (00) - whose amplitudes are 
related by an isospin quadrilateral relation. There are 
nine measurements of these processes that can be made: 
the four branching ratios, the four direct CP asymme- 
tries A CPl and the mixing-induced CP asymmetry S CP in 
B d — > tt°K°. Several years ago, these measurements were 
in disagreement with the standard model (SM), leading 
some authors to posit the existence of a U B — > nK puz- 
zle" [lj]. Subsequently, there was an enormous amount 
of work looking at B — ► irK decays, both within the SM 
and with new physics (NP). 

One simple way to see this discrepancy is to define the 
ratios of branching ratios 
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Within the SM, it is predicted that R c is approximately 
equal to R n 0- QCD-factorization techniques Q can be 
used to take corrections into account. These yield [i| 



f? — 1 lfi+0-22 p 1 ic+unj 
J^n — r.iu_Q , n c — i..ij_Q yj 



However, the pre-ICHEP04 data was Q 
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The values of R n and R c revealed a clear discrepancy 
with the SM. Even the pre-ICHEP06 data was in dis- 
agreement with the SM. 

However, the ICHEP06 data is 
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which is now in agreement with the SM. This has led 
some authors to claim that there is no longer aB-t *kK 
puzzle 0,0]. 

In the present paper, we revisit the question of whether 
or not there is a B — > irK puzzle. As we will see, the B — > 



■nK puzzle is still present, and is in fact more severe than 
before. However, it now is present only in CP-violating 
asymmetries. If one looks at only the measurements of 
the branching ratios, as above, one will not see it. We also 
examine the question of which NP scenarios can account 
for the current B — > ixK data. As such, this paper can 
be considered an update of Ref . 0] , in which the present 
authors were involved. 

We begin by writing the four B — > ttK amplitudes us- 
ing the diagrammatic approach. Within this approach 
Q, the amplitudes can be written in terms of six dia- 
grams: the color-favored and color-suppressed tree am- 
plitudes T 1 and C , the gluonic penguin amplitudes P' tc 
and P' uc , and the color- favored and color-suppressed elec- 
troweak penguin amplitudes P' EW and P' BW - (The primes 



on the amplitudes indicate b 
plitudes are given by 



transitions.) The am- 
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where we have explicitly written the weak-phase depen- 
dence, (including the minus sign from V t * b V ts [-P/J), while 
the diagrams contain strong phases. (The phase infor- 
mation in the Cabibbo-Kobayashi-Maskawa quark mix- 
ing matrix is conventionally parametrized in terms of the 
unitarity triangle, in which the interior (CP-violating) 
angles arc known as a. /J and 7 [9}.) The amplitudes 
for the CP-conjugate processes can be obtained from the 
above by changing the sign of the weak phase (7). 

Within the SM, to a good approximation, the diagrams 
P' BW and P'° w can be related to T" and C using flavor 
SU(3) symmetry [ljj: 
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Here, the Cj are Wilson coefficients and i? = 

(Kt^)/(1^* 6 K S )|. In our fits we take R = 48.9 ± 1.6 

m- 

In Ref . Q , the relative sizes of the B — > 7r A~ diagrams 
were roughly estimated as 



1 : |P/ C | , 0(A) : |T'|, 
0(A 2 ) : \P' B %\ 
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where A ~ 0.2. These estimates are expected to hold 
approximately in the SM. If one ignores the small 0(A 2 ) 
diagrams, the four B — > irK amplitudes take the form 
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Mode 


Bi?[KT 6 ] 


Acp Sep 


B+ 




23.1 ± 1.0 


0.009 ± 0.025 


B+ 


^7T°A+ 


12.8 ±0.6 


0.047 ± 0.026 


R° 
B d 




19.7 ±0.6 


-0.093 ±0.015 


no 


^7r°A° 


10.0 ±0.6 


-0.12 ±0.11 0.33 ±0.21 



TABLE I: Branching ratios, direct CP asymmetries A C p, and 
mixing-induced CP asymmetry Sep (if applicable) for the 
four B — > ttK decay modes. The data is taken from Refs. [7J|] 
and 0]. 



From the above, we see that the amplitudes 
V2A(B+ -> ir°K+) and A(B% -> tt- K+) are equal, up 
to a factor of P' EW . Since \P' EW \ < |P/ C |, we therefore ex- 
pect the direct asymmetry ^4 CP (0±) to be approximately 
equal to A CP ( — h): 
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where <5 T / is the strong-phase difference between T 1 and 
P/ c . However, looking at the ICHEP06 B ^ ttK data in 
Table 1, we see that these asymmetries are very differ- 
ent. Thus, the 2006 B — > 7rAT data cannot be explained 
by the SM. It is only by considering the CP-violating 
asymmetries that one realizes this. 

To investigate this numerically, we have performed 
a fit to the 2006 ICHEP06 B -> ttK data using the 
parametrization of Eq. JS]). We find that Xrnin/d-°-f- = 
25.0/5 (1.4 x 10~ 4 ), indicating an extremely poor fit. 
(The number in parentheses indicates the quality of the 
fit, and depends on Xmin an d d.o.f. individually. 50% or 
more is a good fit; fits which are substantially less than 
50% are poorer.) 



The only way to account for the data within the SM 
is to include C", since it appears in only one of the two 
amplitudes y/2A{B+ -> tt°K+) and A(P° -> n~K+). 
Thus, as a first step to remedy the problem, we follow 
Refs. 15j and add only C [from Eq. |(3J)] to the ampli- 
tudes of Eq. ijHJ). From the approximate formulae 



A CP (0+) 
A CP {-+) 



V 



P' 

r 



p' 



sin S T i sin 7 
sin S T i sin 7 — 2 



p' 



(10) 

sin 5 C > sin 7 , 



where S c / is the strong-phase difference between C and 
P/ c , we see that a large value of \C'\ can give the cor- 
rect sign for A CP ( — h) when sin<5 C ' has a different sign 
from sin^ T '. This is confirmed numerically. A good 
fit is obtained: xL in /d.o.f. = 1.0/3 (80%). We find 
|P'| = 47±1 eV, |T'| = 8.1±3.5eV, \C'\ = 13.0±3.2eV, 
S T , = (154±10)°, 8 C > = (-154±7)°. However, note that 
\C'/T'\ = 1.6±0.3 is required (we stress that correlations 
have been taken into account in obtaining this ratio). 
Since this value is much larger than the naive estimates 
of Eq. ([7]), this shows explicitly that the B — > irK puzzle 
is still present. In Ref. 0] (2004), \C'/T'\ = 1.8 ± 1.0 
was found. We thus see that the puzzle has gotten much 
worse in 2006. (See Ref. [l(| for another approach.) 

We note in passing that the two fits give 7 = (62.5 ± 
11.1)° and 7 = (50.0 ± 5.6)°, respectively. Thus, both 
fits give values for 7 which are consistent with the value 
obtained via a fit to independent measurements: 7 = 



12] 



59.orlf 

We can also perform a fit which incorporates all the 
C(A 2 ) terms in Eq. ([5]) (i.e. P' uc is included). Once again, 
we find a good fit: x* mn / d -°-f- = °- 7 A ( 79 %) ( th c values 
of A CP (0+) and A CP {— ±) are reproduced due to C'-P' uc 
interference). However, a smaller value of \C'\ is found, 
compensated by a large value of | P' uc \ : \C'/T'\ = 0.8±0.1, 
\PiJT'\ = 1.7 ± 0.6. We also find that 7 = (30 ± 7)°. 
Thus, the problem in \C'/T'\ has been alleviated, but we 
encounter new difficulties in |P„ C /P'| and 7. 

Having established that there still is a B — > irK puz- 
zle, we now turn to the question of the type of new 
physics which can be responsible. There are a great 
many NP operators which can contribute to B — > 7rA' 
decays. However, this number can be reduced consid- 
erably. All NP operators in b — > sqq transitions take 
the form O^p ~ sTibqTjq (q = u,d,s,c), where the 
Tij represent Lorentz structures, and color indices are 
suppressed. These operators contribute to the decay 
B — ► ttK through the matrix elements (ttK\ 0$J? \B), 
whose magnitude is taken to be roughly the same size as 
the SM b — * s penguin operators. Each matrix element 
has its own NP weak and strong phase. 

We begin with the model-independent analysis de- 
scribed in Ref. [l7| and used in Ref. • Here it is argued 
that all NP strong phases are negligible. This allows for a 
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great simplification. If one neglects all NP strong phases, 
one can now combine all NP matrix elements into a single 
NP amplitude, with a single weak phase: 



{irK\Of q \B)=A q e 
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Now, B — > ttK decays involve only NP parameters re- 
lated to the quarks u and d. These operators come in two 
classes, differing in their color structure: s a Tib a qpTjqp 
and s a Tibp qpTjq a (q = u,d). The matrix elements of 
these operators can be combined into single NP ampli- 
tudes, denoted A'' q e 1 ^ and _4 /c <V*f , respectively 
Here, §' q and <&'° are the NP weak phases; the strong 
phases are zero. Each of these contributes differently to 
the various B — > ttK decays. In general, A' q ^ A' c ' q and 
<5g 7^ &q ■ Note that, despite the "color-suppressed" in- 
dex C, the matrix elements A lc,q e i ®i are not necessarily 
smaller than the A'^e^t . 

The B — > ttK amplitudes can now be written in terms 
of the SM amplitudes to 0(X) \P' EW and T" are related 
as in Eq. (jBJ)], along with the NP matrix elements (l8j : 



A+ Q = -Pi + A' c d e l ^ C , 
V2A 0+ = Pl c -T'e*-P' BW 
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Note that the above form of NP amplitudes still satis- 
fies the isospin quadrilateral relation. If P' tc is dominant, 
one will also find that R n ~ Rn and the sum rules for CP 
asymmetries given in Ref. [19j are satisfied. Therefore 
the fact that the SM satisfies those sum rules does not 
necessarily rule out this kind of NP. (Of course, if the NP 
is about the same size as P( c , one can find that R„ =/= R c 
and that the sum rules are violatied.) 

The value of 7 is taken from independent measure- 
ments (7 = 59.0^4 g 12|), leaving a total of 10 theoret- 
T\, \A'' comb \, \A'°> U \, \A' c ' d \, 3 

With 



\PL\ 



ical parameters: 

NP weak phases and two relative strong phases 
only 9 experimental measurements, it is not possible to 
perform a fit. It is necessary to make some theoretical 
assumptions. 

As in Ref. Q, we assume that a single NP amplitude 
dominates. We consider the following four possibilities: 
(i) only A'' comh ^ 0, (ii) only A' c ' u ^ 0, (iii) only A'°> d ± 



0, (iv) A 



A' 



A >,comb = o ( isospin _ 



conserving NP). 

However, before presenting the results of the fits, we 
can deduce what the results should yield. Earlier, we saw 
that the SM can reproduce the ICHEP06 B — ► ttK data, 
but only if C is anomalously large. Of the three NP 
matrix elements, only _4'> comb e l * appears only in ampli- 
tudes of decays which receive contributions from C (A 0+ 



reproduce a large C (or P' EW |20j), it is necessary to add 



interference of P' EW 



and A! c,d e 1 ^^ are not necessary. 



This is 



and ,4 00 ); A' a ' u e i9 '" and A IG ' d e W ° appear also in ampli- 
tudes with no C (A +0 and A~ + ). Thus, if one wishes to 

In light of this, we expect a good fit only for possibility 
(i) above; (ii), (iii) and (iv) should have poor fits. (Note 
that fit (ii) might be marginally acceptable due to the 
and A'°' u e w ° .) 
borne out numerically. We find (i) 
xLJd-o.f. = 0.6/3 (90%), (ii) xLJd-o.f. = 
4.0/3 (26%), (iii) xliJd-o.f. = 21.1/3 (0.01%), (iv) 
XmiJd-o.f. = 12.8/4 (1.2%). (Regarding the d.o.f.'s: in 
fits (i)-(iii), we fit to the three magnitudes \P[ C \, \T'\ and 
|^4| (NP), two relative strong phases, and the NP weak 
phase <J>, so the d.o.f. for these three fits is 9 — 6 = 3. 
In fit (iv), one of the NP parameters can be absorbed 
into l-P/J, so the number of d.o.f.'s is increased by one.) 
Thus, as expected, a good fit is found only if the NP is 
in the form of A'' comb e W . 

This is the same conclusion as that found in Ref. 0]. 
Thus, not only is the B — > irK puzzle still present, 
but it is still pointing towatds the same type of NP, 



A''' 
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7^ (this corresponds to NP in the elec- 



troweak penguin amplitude). For this (good) fit, we find 
\T'/P'\ = 0.09, \A'< comb /P'\ = 0.24, $' = 85°. We there- 
fore find that the NP amplitude must be sizeable, with a 
large weak phase. 

If one does not wish to neglect the NP strong phases, 
one can still reduce the number of NP operators. It can 
be shown that an arbitrary amplitude can be written in 
terms of two others with known weak phases. This is 
known as reparametrization invariance (RI) plj . Apply- 
ing this to the NP amplitudes, and taking the two weak 
phases to be and 7, we can write 



{*K\ Of q \B) 
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contain only strong phases. 
ttK amplitudes can now be written as in 
where each of the NP amplitudes is expressed 
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12| , there are 



as in Eq. (JT3J) above. Taking 7 = 59.0 
a total of 15 theoretical parameters: the eight magnitudes 
\PLl in l^;? 6 !' Ka^ Ko^l. and seven relative 
strong phases. With only 9 measurements, we must again 
make theoretical assumptions, and we consider the four 
possibilities (i)-(iv) described above. 

As before, we expect a good fit only for possibility (i). 
For the d.o.f.'s, in the fits there are seven theoretical pa- 
rameters: the four magnitudes \P[ C \, \T'\, \Aq\ and \A q \, 
and three relative strong phases. In fit (ii), ^4 CP (+0) 
and j4 cp (00) are identically zero, and S , CP (00) = sin 2(3. 
Thus, these measurements do not constrain the parame- 
ters; there are only six constraining measurements. This 
is fewer than the number of parameters, so a fit cannot 
be done in this case. The fit can be done in cases (i), (iii), 
and the d.o.f. is 9 — 7 = 2. In fit (iv), Aq can be absorbed 
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into P[ c , so the d.o.f. is 9 — 5 = 4. The results of the fits 
are: (i) X 2 ml Jd.o.f. = 0.5/2 (79%), (iii) xLjd.o.f. = 
18.5/2 (0.01%), (iv) X 2 mi Jd.o.f. = 12.8/4 (1.2%). As 
expected, only fit (i) [NP = A'> comb ] gives a good fit. We 
find \T'/P'\ = 0.09, \Al/P'\ = 0.15, \A V /P'\ = 0.28. We 
therefore find that the NP amplitude must be sizeable. 

We have considered two different parametrizations of 
the NP effects: negligible NP strong phases, and the gen- 
eral case. Although the results of the fits are not numer- 
ically identical, the conclusions are similar. 

Of the four new-physics models examined in this paper, 
only one produces a good fit to the 2006 B — > ttK data: 
case (i), y\' ,comb ^ 0. The quality of fit is extremely good 
here. The NP amplitude must be sizeable, with a large 
weak phase. Fit (ii) (A' c ' u ^= 0) is marginal if the NP 
strong phases are negligible (the fit cannot be done in 
the general case). Here A CP (0+) « A CP { — h) is found 
due to the interference of P' EW and *4' c,t V*" . Fits (iii) 
(A' c,d ^ 0) and (iv) (isospin-conserving NP) yield very 
poor fits, and are ruled out. 

As mentioned in Ref . Q , we note that we have assumed 
that one NP operator dominates. However, any specific 
NP model will generally have more than one effective NP 
operator, and so the more general case might be used to 
explain the B — > irK data. 

In conclusion, in this paper we have done a fit to the 
2006 B — * irK data to see if new physics (NP) was still 
indicated. Note that a fit to the full data is necessary; 
if a subset of the measurements is used, erroneous con- 
clusions can be reached. If only the "large" diagrams of 
the SM are included in the amplitudes, a very poor fit is 
found. If the "small" diagram C is added, a good fit is 
obtained. However, \C /T'\ = 1.6 ±0.3 is required. The 
fact that this ratio is large indicates that the B — ► irK 
puzzle is still present. The fact that the error is so small 
shows that the puzzle has gotten much worse in 2006. 

We have also examined which type of NP can be re- 
sponsible for the B — ► irK puzzle. We assume that one 
NP operator dominates, and consider (i) only A'' com ^ 
0, (ii) only A' c ' u ^ 0, (iii) only A' c > d ^ 0, (iv) isospin- 
conserving NP: A' c ' u e i9 '" = A lc > d e w i , A'> comb = 0. Of 
these, the only (very) good fit is found for model (i). It 
corresponds to a modification of the electroweak penguin 
amplitude. Note that these results are the same as those 
found in Ref. Q. This shows that, although the puzzle 
is worse than in 2004, its resolution is the same. 
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